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a 

Speed of sound 

R 

Gas constant (for nitrogen 


Constants of the Beattle- 


- 896.82 Nm/kg K) 

1 

Dridgeman equation 
(1 - 1,2. . .5) 

s 

entropy 

^A 

Lift coefficient 

T 

Temperature 

n 

c_ 

Pressure coefficient 

t 

Thickness of pi’oflle 

p 

(* critical pressure co- 


Velocity of the undisturbed 


efficient) 


oncoming flow 

Cp 

Specific heat at constant 

u,v 

Disturbance velocities 

p 

pressure 

x,y 

Coordinates 

c 

Specific heat at constant 


Angle of attack, expansion 

V 

volume 

a 

coefficient 

1 

Enthalpy 

r 

Hayes' similarity parameter 

1 

Depth of profile 


Similarity parameter 

Ma_ 

00 

Mach number of undisturbed 
oncoming flow 

< 

Ratio of specific heats 

f'l 

Constants of the Jacobsen 
equation (1 = 1,2... 5) 

p 

P V 
Density 



<]> , ^ 

Potential 

p 

pressure 
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1. Introduction 

The Reynolds number is an important simulation parameter 

for applying results obtained from small-scale wind tunnel models 

to the full scale model. Since especially in the transonic 

speed range the laws of similarity are not sufficiently known, 

very serious efforst have been made for years to Increase the 

Reynolds number in wind tunnel experiments. The two traditional 

ways of increasing the Reynolds number — namely increasing 

the static pressure or increasing the linear dimensions — soon 

led to limits determined by costs and the sturdiness of the model. 

Very early on it was recognized by R. Smelt [Ij that a reduction 

in temperature could bring about a considerable improvement. 

This idea was first taken up again in 1972 by M. J. Goodyer and 

R. A. Kilgore [2,33* In a pilot tunnel constructed by 

NASA at Lnagley Field the technical feasibility of wind tunnel 

experiments at rest temperatures of 90K-120K was impressively 

demonstrated. However, since the flow medium (as a rule nitrogen) 

no longer behaves like an ideal gas in the vicinity of the 

llquifactlon limit, the applicability of test values again came 

into question. Studies by J. Adcock, R. A. Kilgore and 

E. J. Ray [4] showed, however, that in the generally used gas- 

dynamic calculations the deviations from the behavior of an 

ideal gas are only very small. In particular, from values 

tabulated by R. T. Jacobsen and R. B. Stewart [53, Adcock 

2 

determined a pseudo adiabatic exponent, a « a p/p, for typical 
experimental coridltlons. This deviated only slightly from a 
value of 1.40. Further studies, especially by H. Vivland [6J 
soon showed that neither the ratio of the specific heats tc nor 
the pseudo adiabatic exponent a is a determining factor for the 
similarity of the flow. The following value was given by Vlviand: 

K = 1 + (da“/di) 

s 

this replaces the < in the transonic potential equation. This 
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essentially universal similarity parameter was given by 

W. D. Hayes [7] as early as 195^, however In a somewhat different 

f 0 rm : 

(F » (x* + l)/2 ■ (dpa/dp) */a) . 

Studies on the parameter f, especially with respect to compressible 
flows of real gases and liquids, were carried out by 
P. A. Thompson [8], P. A. Thompson and K. C. Lambrakls [9] and 
P. A. Thompson and B. A. Sullivan [10]. 

In this paper, after a brief refresher discussion of how 
<* Is derived [6], an equation will be giver which permits a 
simple calculation of k* from thermal and caloric parameters of 
state. Numerical values for are obtained by means of the 
Beattle-Brldgeman equation of state [11] and from [5] for 
nitrogen as the working medium. Applications of the Beattle- 
Brldgeman equation to various transonic flows are given by 
B. W. Wagner and W. Schmidt [12]. The influence of the. similarity 
parameter on transonic flows Is shown here by the flow around a 
supercritical profile. 

2. Derivation of /k 

Several ways of deriving tc* will be discussed In this section 
section. In so doing, we will follow the line of thought of 
VI viand [6]. 


( 1 ) 


( 2 ) 

OUIGLNAL PAGK l3 
OF POOR QUALITY 


The following initial equations are used: 
the potential equation 


the energy equation 


i f = i 4 /2 

' X V ' cr* co' 



and an equation of atate 


2 2 , . 

a ■ a (i, s) 


( 3 ) 


Using the theory of small disturbances 


4' "U +u , ❖ • V 

X CO y 


a ~ a 


UJ 
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one obtains the following equation for Isentropic flows: 

2 2 ,^ 2 ,,. 


oo 


( 5 ) 


and with Eq. (2) one obtains the following equation 

2 .2. 


( ht ) ■■-(%-) 


2 u_ 

CD 


( 6 ) 


Applying the theory of small disturbances to the first 
factor of Eq. (1) gives us the following: 


( 7 ) 


ai 


CO 


Substituting Eq. (6), we get 
^ 2 


.( .Ji) « 1 . - Ma^ ^2 + Ma“ (x* - ^))r^ (8) 

\ a / to co\^ CO / 


with 


* 

X = 1 ♦ 


(^) 


( 9 ) ^ 


Thus the potential equation reads as follows: 

f 1 - Ma^ - Ma^ (2 » Mn^ (x* - D)oJP^.^. ^ O,.,. = 0 
\ CO CO' CO • X/ x\ >y 


( 10 ) 
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This equation lo formally Identical to the transonic 
potential equation for Ideal gases. In the case of the real 
gas only the ratio of the specific heat must be replaced by 
K« according to Eq. (9). 

In dealing with tabular figures (e.g. [5])i the num*» *lcal 
determination of <* runs Into considerable dlffl . .1 ties . There- 
fore, In the following section the attempt will be made wO find 
a more convlenent analytical form for 


The first law of thermodynamics yields the following for 
Isentroplc changes of state: 


pdi ‘ dp , 

so that Eq. (9) can be written In the following form: 
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If one takes a as a function of p and T, one obtains 


( 11 ) 


(12) 




with the Identity (of. E. Schmidt [13]» ch. 11) 



one obtains for Eq. (13) 




(13) 


( 1 ^) 


(15) 


(16) /6 



For an ideal gas one Inunedlately obtains 



p 

Assuming that a Is a function of the thermal parameters of 
state p and T, another form for can be determined In a 
similar manner which has been given by E. Wedemeyer [14], This 
reads as follows: 



Eq. (16) was developed In order to be able to determine *«• 
from tables In which the calculated parameters of state are 
given as a function of p and T. Eq. (18), by contrast, Is 
directed more towards calculation from equations of state of the 
form p ■ f(p,T). In the next section a few ways of calculating 
K* from equations of state will be demonstrated. The numerical 
values determined from tables [5] are also given. 


3. Calculation of k* 

3.1 General Remarks 

For the most part, equations of state are given In the form 

p « pip, T) 

(19) 

In order to calculate k* It is necessary to know the speed of 
sound a , the specific heats Cp and c^ and the corresponding 
derivatives of the pressure and the speed of sound. To this 
end the following equations are used [5,13]: 
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c 

V 


(3|./JT)^ 

7 


2 

a 





( 20 ) 


( 21 ) 

( 22 ) 


3.2 The Beattie- Bridtj:eman Equation of State 

The Deattle-Bridgemean equatlor, of state wan given for the 
first time in 1927 [11]. It can be used for many gases. Up to 
gas densities of P ■ 1^0 kg/m^ the gas conditions calculated 
with this equation are in excellent agreement with observations. 
It has 5 constants and is written as follows: 


\t - p-H*T(l - pA^/r^)(l A^p(l - A^o)) - Ajp^(l - A^p) . (23) 

The constants are taken from [11]. In the MKS system for 
nitrogen they are as follows: 

A, « 173. G A„ * 1.8013.10’^ 

1 J 

.4 -4 

A,, -^9, 342-10 , A . » - 2.4U7-10 . A. * MOy.o 

*5 • 4 3 

(P(ki;/m^l , K(Nm/kgK| . T(KI , pfN/m^J) . 

The necessary derivatives can be given immediately: 

(!t) " P-»(1 +2pAj/T^)(l tAjPd . A,_P)) (2^) 


1 . 
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2* 


II = 7-11.99 Km/kt» K 


( for nitrogen) . 


(1^) - U-T(l - p/\j^(2 f - 3pA^ - 4p^A^A^/T^)) (2‘j) 


- AjP( 2 - 3A^p) 

\ ■<‘-v>id*“'”'V‘ -A2P(i-ipA,))/T’ 


• (^t) 


3.^ The Equation of State after R. T. Jacobsen /fl 

The tables by Jacobsen [5] are based on an equation of 
state with 39 constants which were determined from values 
measured in a pressure ranpe of p ■ 0.1 to 10,000 bar and a 
temperature range of T ■ 63 to 2,000 K. However, for the 
operation of a cryogenic wind tunnel the olanned pressure range 
and temperature range are much smaller (p up to 8 bur, T up to 
300 K) so that here It will be attempted to get by only with 
5 constants. Such an equation, accoi'dlng to Iti], runs as 
follows : 


pKT -» + N^\"T + N., ♦ N.^/T ^ 


(27) 


In the HK3 system the constants for nitrogen are as follows: 


Nj ' 0. 1V572 

N , -0. -HOHO-IO 
4 

(p(kii/tn^l 


N « 0. 1381}}.in‘ 

» - 0. b4(iH-ll/' , 

D 

HlNin/ktiKI . T(K) 



n. 31 49. '■.•10^ 




The simple derivatives turn out to be the following: 




( 28 ) 
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(29) 


(fp) ‘ + N^/T ^ Nj^/T^) 


c , • (t\ ) + P(0. 2r> N JVT’- 2iN’ /T^ - 6N,/T^) . 


id 


(3D 


2 

3D Calculation of < [sic] 

In order to calculate values for <* from the equations of 
state given in section 3*2 and 3*3 the following procedure is 
selected. For a given pair of values (p, T) the density p is 
determined Iteratively from Eq. (23) or (27). Then these values 
are determined: 


2 2 

The derivatives <^a /fit and fia /fip are determined numerically 
with the help of Spline polynomials. Then can be calculated 
directly from Eq. (16). Some results are shown in the figures 
given at the end of this paper. Figs. 1 and 2 show the derivatives 
of the speed of sound according to temperature and pressure based 
on equations of state (23) and (27). If one formed the cor- 
responding derivatlv'js from [5]* then deviations appear which 
suggest that the table is not sufficiently precise. Therefore, 

It is too Inaccurate for calculating <*. The perceptible dif- 
ferences between Beattle-Bridgeman (Eq. (23)) and Jacobsen 
(Eq. ^27)) clearly show that even these equations do not yet 
describe the gas conditions with sufficient accuracy. 

rTg. 3 shows the values for k* calculated from the two 
equations of state In comparison with a few values determined 
directly from the table. 
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The following conclusions can be drawn from Fig. 3: 

- In the Interesting range for a cryogenic wind tunnel the 
similarity parameter »c* deviates only slightly from 1.^0. 

- At T » 105 K nitrogen behaves like an Ideal gas. 

- The real gas effects are more strongly Influenced by the 
pressure than by the temperature. 

- Differences depending on the equation of state selected 
are manifested only in the vicinity of the saturation 
point . 

- The values determined from the table are not accurate 
enough. 

Applications 

In conclusion, we will briefly discuss the possible Im- 
plications of a change in the similarity parameter. 

With the example chosen here, the pressure distribution 
with a mach number Ma„ « 0.75 for the oncoming flow and an 
angle of attack of a ■ 0® Is calculated according to the method 
of Garabedlan and Korn [16,17] for a supercritical profile 
which was developed by Messerschmltt-Bfllkow-Blohm and Co. [15l* 

In so doing, the similarity parameter k* Is varied within the 
limits encountered In practice (Fig. 3)* Fig. 4 shows the 
pressure distribution for • 1.4. Fig. 5 shows the 
recorapression range in which the k* Influence is detectable. It 
is obvious that the effect of a variation In <* is small both 
with respect to the pressure distribution and with respect to 
lift. 

!; . Summary 

A formula Is given for the determining parameter k* « 1 + 

P 

(6a /6i)_ used for estimating real gas effects. The formula makes 
It possible to determine this parameter both from tables and 
from equations of state (p(p,T)). The equa’-ions of state of 


/lO 
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Beattle-Brldgeman and Jacobsen are applied to the formula. Neither 
equation, however, permits a sufficiently exact description of 
the gas conditions. Using tne calculated pressure distribution 
of a supercritical profile, the effect of a variation In k * 
within the limits encountered In practice Is demonstrated. 

The results obtained permit the conclusion that nitrogen. In 
the pressure range considered and In the vicinity of the 
saturation point shows only slight real gas effects. 
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Fig. 1. 6a /6T at constant 
pressure (for nitrogen) 

Key: A) According to equation 



temperature (for nitrogen) 
Key: A) According to equation 



Fig. 3. * as a function of 

temperature and pressure (for 
nitrogen) . 

Key: A) Equation 

B) From tabulated values 

C) Saturation 
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Fig. 4. Calculated pressure distribution 
on a supercritical profile MBB-A3 with 
K* =» 1.40 (ideal gas) 

Key: A) Zone of influence, see Fig. 5 
B) Profile 





Fig. 5. Pressure distribution on the 
suction side of the r*3B-A3 profile 
(segment) for various values. 



